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We present a scaling theory for the entanglement spectrum under an external driving. Based on 
the static scaling of the Schmidt gap and the theory of finite-time scaling, we show that the Schmidt 
gap can signal the critical point and be used to estimate the critical exponents no matter in the 
hnite-size scaling region or in the finite-time scaling region. Crossover between the two regions is 
also demonstrated. We verify our theory using both the one-dimensional transverse-field Ising model 
and the one-dimensional quantum Potts model. Our results confirm that the Schmidt gap can be 
regarded as a supplement to the local order parameter. 
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I. INTRODUCTION 

Over the past decades, studies on the properties of 
entanglement have provided deep insight in understand¬ 
ing quantum matters and quantum phase transitions [1- 
5]. The entanglement is usually measured by an entan¬ 
glement entropy and an entanglement spectrum [6-10]. 
Both the entanglement entropy and the entanglement 
spectrum can be used to characterize the phase of quan¬ 
tum matters. For example, a topological entanglement 
entropy is proposed to identify the intrinsic topologi¬ 
cal order [11-13], while the degeneracy of the entangle¬ 
ment spectrum is used to classify the symmetry-protected 
topological order in one dimensional (ID) spin chains 
[14, 15]. Aside from the applications in characterizing 
the phases of matter, the scaling behavior of the entan¬ 
glement [16-20] near a quantum critical point is proposed 
to locate the quantum critical point and estimate the crit¬ 
ical exponents [21]. As a significant characterization in 
condensed matter physics, entanglement can be detected 
in experiments and numerical simulations [22]. 

On the other hand, non-equilibrium dynamics of quan¬ 
tum many-body systems has also attracted attention 
[23, 24]. These studies are partly motivated by the 
progress of technologies in detecting and manipulating 
dynamic phenomena in cold atom experiments [25, 26]. 
For the driving dynamics (quench) of changing the dis¬ 
tance g to a critical point with a rate R, the Kibble- 
Zurek mechanism (KZM) divides the evolution into 
two adiabatic stages and an intervening impulse stage 
[23, 24, 27, 28]. In the impulse stage, according to the 
theory of finite-time scaling (FTS) [29-32], the external 
time scale C,d which is induced by the exter¬ 

nal driving, dominates the dynamics, since it is smaller 
than the intrinsic time scale C,s ~ where v and 2 

are critical exponents and r = z + l/iz. The FTS theory 
further reveals that macroscopic quantities are related to 
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the driving rate R with some exponents. For example, 
the number of topological defects n, which are produced 
as a result of the diabatic process in the impulse stage, 
is scaled as n ^ 7?“^/’'. This is consistent with the pre¬ 
diction of the KZM [27, 28]. 

As the entanglement provides another way to under¬ 
stand quantum criticality besides local operators, it is 
meaningful to investigate its nonequilibrium dynamic 
scaling. In the driving dynamics, it has been shown 
that the entanglement entropy is proportion to the loga¬ 
rithm of the driving rate at the critical point [23, 24, 33]. 
In real-time relaxation dynamics, the entanglement en¬ 
tropy increases linearly with time at the short-time stage 
[34, 35], while in imaginary-time relaxation, it increases 
as a logarithm of time at the critical point [36]. Besides 
investigations on the entanglement entropy, the dynamic 
scaling of the entanglement spectrum has also been stud¬ 
ied in both driving and relaxation dynamics [37, 38]. In 
the transverse-field Ising model, for the driving dynam¬ 
ics on which we shall focus in the following, it has been 
shown that at the critical point, the Schmidt gap, AA, de¬ 
fined as the difference between the two largest entangle¬ 
ment spectra, increases with the driving rate as a power 
law AA ~ 77“, with a being estimated to be 0.064 [37]. 
In addition, for large driving rates, no finite size effects 
have been observed [37] . 

Questions arisen are then: (i) How to characterize the 
scaling behavior of the Schmidt gap in the whole driv¬ 
ing process? (ii) Whether is the exponent a, reported in 
Ref. 37, a new exponent or a combination of other known 
exponents? (iii) Is there a unified scaling theory to char¬ 
acterize the driving dynamics for large and small driving 
rates? (iv) Are these results applicable in other models? 

To answer these questions, in this paper, we propose 
a dynamic scaling theory to describe the universal dy¬ 
namics of the Schmidt gap according to the theory of 
FTS. According to the theory, we shall show that the ex¬ 
ponent a is not an independent exponent but rather is 
Plvr, where /3 is the critical exponent of order param¬ 
eter. The crossover behavior between the FTS region 
and the hnite-size scaling (FSS) region is also studied. 
A dynamic method to determine the critical point and 
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critical exponents is proposed. To confirm the validity 
and show the universality of the scaling theory, both the 
ID transverse-field Ising model and the ID three-state 
quantum Potts model are employed. 

The rest of the paper is organized as follows. After 
introducing briefly the Schmidt gap and its static scaling 
behavior in Sec. II, we present the dynamic scaling theory 
of the Schmidt gap under an external driving according to 
the FTS theory in Sec. III. Then, in Sec. IV we introduce 
the models and our methods, while in Sec. V, we verify 
the theory by numerical simulations of the driving critical 
dynamics for the models. Finally, a summary is given in 
Sec. VI. 


II. BRIEF REVIEW OF THE SCHMIDT GAP 
AND ITS STATIC SCALING 

A. Definition of the Schmidt gap 

We first introduce the definition of the Schmidt gap. 
An arbitrary wave function |'0) can be represented into 
a bipartite form by the Schmidt decomposition as 

I^A) = ® (1) 

k 

where \ijj\) )) is the fcth eigenvector of the reduced 
density matrix of the left (right) subsystem and Afc, the 
fcth eigenvalue of the reduced density matrix, represents 
the entanglement spectrum and satisfies Afc > 0 and 
J2k^k — 1- IV'fe) and IV'fc') are orthogonal to each other 
and Afc is identical for the reduced density matrix of the 
left and right subsystems. For simplicity, we assume the 
array of Afc is arranged in a descending order with k. As 
a result, the Schmidt gap, AA, is defined as [18, 19] 

AA = Ai - A2. (2) 

For an infinite system, which exhibits translational sym¬ 
metry, AA does not depend on the location of the lattice, 
while for a finite system, AA does. Accordingly, for a fi¬ 
nite system, we choose AA as the one at the center point 
of the system. This choice can weaken the boundary ef¬ 
fects. 


B. Static properties of the Schmidt gap in 
quantum phase transitions 

It has been proposed that the Schmidt gap can be re¬ 
garded as an order parameter to characterize quantum 
phases and quantum phase transitions [18, 19]. To fur¬ 
ther confirm it, we compare the behavior of AA with the 
local order parameter, M, defined, for example, as the lo¬ 
cal magnetization in the transverse-field Ising model [41]. 
It has been shown that in the disordered phase, M = 0 
but AA ^ 0, while in the ordered phase, M ^ 0 but 


A A = 0. This means that A A is a complementary quan¬ 
tity to the local order parameter to describe the different 
phases [18, 19]. 

The static scaling behavior of AA can be obtained as 
follows. Note that, it has been shown that the critical 
exponent related to M is identical with that to AA [18, 
19]. In other words, under a scale transformation with a 
rescaling factor 6, AA changes to AXb ^^'^, which is similar 
to the scale transformation of M. Accordingly, near the 
critical point, the full scale transformation of AA for a 
system of size L with a distance g to its critical point 
reads as [18, 19] 

AA(g, L) = b-^/''A\igb^/'',Lb-^). (3) 

Different scaling forms have been obtained from 
Eq. (3). When \g\~’^ L, the system cannot “feel” 
its size and the finite g dominates the scaling behavior. 
In this situation, by setting gb^^'^ = 1, one obtains the 
scaling form [18, 19] 

AX{g,L)=g^f,iLg-'), (4) 

where fi is the scaling function with i being an integer. 
Conversely, when L <g; the correlation length is 

truncated by the size of the system and L dominates the 
scaling behavior. The FSS form, 

AX{g,L) = L-f^/-'h{gL^/''), (5) 

then follows in a similar way [18, 19]. Equations (4) and 
(5) have been verified in the ID transverse-field Ising 
model [18, 19]. 

HI. DYNAMIC SCALING THEORY OF THE 
SCHMIDT GAP IN DRIVING CRITICAL 
DYNAMICS 

In this section, we develop the scaling theory describ¬ 
ing the scaling behavior of the Schmidt gap in driving 
dynamics according to the theory of FTS [29, 30]. For 
g = Rt, the competition among the driving time scale (^d 
and the other two time scales ~ 5”"^^ and Cl 
induced by finite g and L, respectively, controls the driv¬ 
ing critical dynamics near the critical point. Including R 
in Eq. (3), we postulate the scale transformation of the 
Schmidt gap as [29] 

AA(g, A, R) = b-^!''AX{gb'^l'', Afe”). (6) 

Then, similar to the discussion in Sec. IIB, different 
scaling forms in different regions can be obtained by 
choosing suitable scale factors. By setting RW = 1, one 
obtains the scaling form of FTS [29, 30], 

AA(g, A, R) = A-IA-A”). (7) 

In this region, the system falls out of equilibrium and its 
evolution is dominated by the external driving rate, while 
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the finite-size effects are only perturbation. In particular, 
in the thermodynamic limit L oo, we obtain 

The FTS form of the Schmidt gap can also applied to 
determine the critical point and critical exponents similar 
to the order parameter [29-31]. Indeed, from Eq. (8), one 
readily finds that when AA(5, R) = 0, the scaled variable 
is a constant. To be explicit, for the transverse- 
field Ising model, g = — h^c and the transverse field 

values, /ixo, at which AA = 0, satisfy 

hxo = ^xc + (9) 

where c is a constant. Finding for various R and 
fitting the results to Eq. (9), one can estimate the critical 
point hxc and the critical exponent l/vr. Then at the 
critical point, from Eq. (8), A A oc Accordingly, 

we can determine the critical exponent fijvr. 

Similarly, by setting Lb~^ = 1, Eq. (6) leads to the 
FSS form [39, 40], 

AA(5, L, R) = L-P/''h{gL^/'^, RL^), (10) 

which returns to Eq. (5) for R = 0. 

Note that both /a and /s can describe the scaling be¬ 
havior in the FTS and the FSS region. As a result, /a 
and /a fulfill 

f^{x, Y) = (ii) 

or conversely, 

/a(A, Y) = Y-^). (12) 

So, at the critical point, the crossover between the 
FTS region and FSS region can be identical. When 
L~^R~^/'^ <C 1, the system is in the FTS region, and 
A A behaves as [29, 30, 39, 40] 

AA(T, R) = R^/'''^h{L-^R-^/''), (13) 

according to Eq. (7) where f^iY) = /^{O.Y). Equa¬ 
tion (13) shows that AX{L,R) oc when L~^R~^/’^ 

can be neglected. Comparing the scaling behavior in 
Ref. 37, in which AA is proportional to for large sizes, 
we find thus a = jijvr. In contrast, when LTR <g; 1, the 
system is in the FSS region, and AA satisfies 

AX{L,R) = L-^/''friL-R), (14) 

and AA oc in this region when UR can be ne¬ 

glected, where /7(E) = /5(0,E). Moreover, according to 
Eq. (11), in the FSS form, the FTS region behaves dis¬ 
tinctly from the FSS region, because asymptotically, /y 
tends to [40] 

friUR) - (15) 

On the other hand, in the FTS form, the FSS region 
behaves asymptotically as 

feiL-^R-^^’') - (16) 

according to Eq. (12) [40]. These results will be born out 
by simulations shown in Figs. 3 and 4 for the Ising model 
and Fig. 7 for the Potts model below. 


IV. MODEL AND NUMERICAL METHOD 

The Hamiltonian of the transverse-field Ising model is 
[41] 

L-1 L 

= (17) 

i—1 2—1 

where erf and erf are the Pauli matrix in the x and z 
directions, respectively, at site i. The critical point of 
model (17) is hx = hxc = 1 [41]. The exact critical 
exponents are /3 = 1/8, v = 1, and z = 1 [41]. This 
model has been realized in CoNb206 [42]. 

In order to confirm the universality of the scaling prop¬ 
erties of the entanglement spectrum, we also employ the 
three-state quantum Potts model [43] with the Hamilto¬ 
nian 

Hp = - xi E - /lx E ( 18 ) 

2—1 q—1 i—1 

where and are the three Potts spin matrices 


'2 0 o' 


'O 1 o' 

0-10 

= 

0 0 1 

0 0-1 


1 0 0 


and = (M^’^)^. The critical point of model (18) is 

again = hxc = 1, and the conjectured critical expo¬ 
nents are /? = 1/9, = 5/6, and z = 1 [43]. 

We use the time-evolving block decimation (TEBD) 
algorithm [44] to calculate the evolution of the Schmidt 
gap under an external driving. After expanding the wave 
function into a matrix product form via the Vidal de¬ 
composition [45] , each site is attached by a matrix. This 
matrix is updated by acting the local evolution operator, 
which is the Suzuki-Trotter decomposition of exp{—i'Ht). 
By taking the translational symmetry into account, a 
variant of the TEBD algotithm has been developed to 
calculate the dynamics of a quantum chain with infinite 
length. This method is called infinite time-evolving block 
decimation (ITEBD) algorithm [46] . We choose the time 
interval to be 0.01 and the number of the kept states to 
be 50. Three decimal places are kept in our results. More 
accurate results are expected by choosing larger numbers 
of truncation and smaller time intervals. 

We start with the ground state of the Hamiltonian for a 
sufficiently large g and decrease it linearly with time with 
a given R. This is just the reverse of the parallel study 
of the order parameter, in which g is increased linearly 
[31], in accordance with their complementary character. 
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FIG. 1: (Color online) Curves of A A versus R collapse onto 
each other for model (17). The corresponding curves before 
rescaling are plotted in the inset. 


V. VERIFICATION OF THE DYNAMIC 
SCALING 

A. Dynamic scaling of the Schmidt gap in the 
quantum Ising model 

In this section, we shall verify the FTS forms of the 
Schmidt gap for the ID transverse-field Ising model. We 
firstly verify the full scaling form of FTS in Eq. (8) for 
the infinite size situation. Figure 1 shows the evolutions 
of the Schmidt gap for various driving rates. Near the 
critical point, the curves collapse onto each other after 
rescaling with the exact critical point and critical expo¬ 
nents as input. This result confirms the FTS form of 
Eq. (8) and shows that AA has the same scaling dimen¬ 
sion as the local order parameter M even in the nonequi¬ 
librium situation. 

Eigure 2 shows the procedure to determine the critical 
properties according to Eq. (9). The critical point h^c is 
estimated to be h^c = 1.000, in good agreement with the 
exact value, with the critical exponents Xjvr = 0.506 and 
jjjvr = 0.060. By substituting z = 1, these exponents 
lead to /3 = 0.118 and ly = 0.976, which are close to the 
exact values /3 = 1 / 8 and ly = 1. 

Then we focus on the the finite-size effects and verify 
Eq. (13). At the critical point, the dependence of AA on 
for various R is shown in Fig. 3(a). From Fig. 3(a), 
we see that for large L, the curves are independent of L. 
In contrast, for small L, the curves for different driving 
rates collapse onto each other and are independent of 
the driving rate. This is because in the former case L 
is bigger than and FTS dominates, whereas in the 
latter case, their relation is reversed, and the evolution 
is dominated by the finite-size effects. 

Figure 3(b) shows that the rescaled curves collapse 
onto each other, confirming Eq. (13). From Fig. 3(b), 
we find that when is small, the rescaled value. 


FIG. 2: (Color online) Estimation of the critical point and 
critical exponents for model (17). (a) Fitting of /ixo for the 
critical point /ixc and the critical exponent 1/i/r according to 
Eq. (9). (b) Fitting of AA(|(;| = 0) for the critical exponent 
P jvr according to Eq. (8). 


'N 
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FIG. 3: (Color online) (a) AA versus L ^ for various R at the 
critical point of model (17). (b) AAi?“'^A^ versus L~^Rr^l''^ 
according to Eq. (13). The two different regions of FSS and 
FTS are marked. Lines connecting symbols are only a guide 
to the eye. Double logarithmic scales are used. 


AAi?”^/'^’’ is almost a constant. This constant is just 
/5(0) in Eq. (13). Eor large the system enters 

the FSS region, in which the curve is almost a straight 
line with a slope of 0.124 in the double-logarithm scales. 
This confirms Eq. (16), which predicts the power-law re¬ 
lation in the FSS region with an exponent ol^jv [39, 40] , 
whose exact value is 0.125. 

We can also study the crossover from the FSS form 
according to Eq. (14). Figure 4(a) shows that for small 
L, the system is in the FSS region, in which AA is al¬ 
most a constant; while for large L, the system is in the 
FTS region and the curves collapse onto each other again. 
Figure 4(b) shows that the rescaled curves collapse onto 
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FIG. 4: (Color online) (a) AA versus R for various L at the 
critical point of model (17). (b) versus RL^ according 

to Eq. (14). The two different regions of FSS and FTS are 
marked. Lines connecting symbols are only a guide to the 
eye. Double logarithmic scales are used. 
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FIG. 5: (Color online) AA versus R collapse onto each other 
for model (18). The corresponding curves before rescaling are 
plotted in the inset. 


each other and confirms Eq. (14). In the FTS region, 
power function fitting gives the exponent to be 0.062, 
which is close to the exact value Pjvr = 0.0625. This 
result confirms Eq. (15). Note that this is also the slope 
of the asymptotic straight line for large L in Fig. 4(a) 
and is thus a in Ref. 37. In fact. Fig. 4(a) is just Fig. 8 
of Ref. 37, albeit with smaller L and R. Therefore, our 
scaling theory explains well the results of Ref. 37. 


B. Dynamic scaling of the Schmidt gap in the 
quantum Potts model 

To further confirm the scaling theory, we study the 
driving dynamics of the quantum Potts model (18). 
First, we consider the infinite-size situation. Figure 5 
shows that the curves for different driving rates collapse 
onto each other after rescaling by substituting the critical 
exponents for the quantum Potts model. This result con¬ 
firms Eq. ( 8 ) and demonstrates the universality of Eq. ( 8 ) 
in the driving dynamics. Comparing Fig. 5 with Fig. 1, 
one finds that in Fig. 1, the Schmidt gap has some dis¬ 
crete points in which AA = 0, while in Fig. 5, there is a 
continuous region corresponding to AA = 0. The reason 
for this difference may be the symmetry group for the 
quantum Potts model is Z 3 while for the Ising model is 
Z 2 . In spite of this difference, the scaling theory can well 
describe the driving dynamics for both models. 

Similar to the situation in the quantum Ising model, 
we can also use the method according to Eq. ( 8 ) to deter¬ 
mine the critical point and the critical exponents. The 
only difference is that we should choose the points just 
at which AA becomes zero as the “characteristic” points 
to determine the critical point, since the zeros for the 
Potts model are continuously distributed. Figure 6 gives 



FIG. 6: (Golor online) Estimation of the critical point and 
critical exponents for model (18). (a) Fitting of /ixo for the 
critical point ftxc and the critical exponent 1/vr according to 
Eq. (9). (b) Fitting of AA(|g| = 0) for the critical exponent 
j3/vr according to Eq. (8). 


the fitting results. The critical point is estimated to be 
hxc = 1.000 with the critical exponents Xjvz = 0.546 and 
/3/i/r = 0.0576. These exponents result in j3 = 0.106 and 
V = 0.832 by substituting z = 1. They are all close to 
their exact and conjectured values. 

Next, we study the finite-size effects in the driving dy¬ 
namics near the critical point. Figure 7(al) shows the de¬ 
pendence of A A on i? for various L. The rescaling curves, 
obtained by rescaling the data with the size L, collapse 
onto each other well, as displayed in Fig. 7(a2). Similar 
to the situation in the transverse-field Ising model, for 
small RL”, the system is dominated by the finite-size ef¬ 
fects, and AAL^/'^ tends to a constant, indicating that 
AA oc . For large i?L”, the system is in the FTS 

region. According to Eq. (15), the slope is Pjvr. This 
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FIG. 7; (Color online) FSS and FTS and their crossover of 
the Schmidt gap for the Potts model (18). (al) Dependence of 
AA on R for various L given in (a2). (a2) The rescaled curves 
from (al) according to Eq. (14). (bl) Dependence of AA on 
L~^ for various R given in (b2). (b2) The rescaled data from 
(bl) according to Eq. (13). Lines connecting symbols are only 
a guide to the eye. 


exponent is fitted to be 0.060, which is close to the exact 
value 0.061. 


We can also present the above results in the FTS form. 
Figure 7(bl) shows the dependence of AA on L for various 
R, while Fig. 7(b2) displays the collapses after rescaling. 
The slope of the FSS region is estimated to be 0.112, 
which is close to the conjected value of filv. These results 
confirm the scaling analysis in Eq. (16). 


VI. SUMMARY 

We have studied the scaling behavior of the entangle¬ 
ment spectrum in the driving dynamics. A scaling the¬ 
ory based on finite-time scaling has been developed to 
describe the scaling behavior of the Schmidt gap of a 
system under an external driving. We have found that 
the Schmidt gap develops similar scaling behavior as the 
local order parameter does. This shows that the Schmidt 
gap can be regarded as a supplement to the local order 
parameter. To illustrate it, we have utilized the FTS 
form of the Schmidt gap to determine the critical point 
and critical exponents. The crossover behavior between 
the FSS region and the FTS region has also been explored 
and the previously published results have been well ex¬ 
plained within the present theory. The universality of the 
scaling theory are confirmed in both the transverse-field 
Ising model and the three-state Potts chain. 
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